Epidemic Dynamics of Interacting Two Particle Species on Scale-free Networks 



O 
O 

(N 

wo: 
< 

(N 



o 



CO 



I 

C 

O 

o 



oo 
o 

O 

c3 



o 
o 



X 



Yong-Yeol Ahn, 1 Naoki Masuda, 2 Hawoong Jeong, 1 and Jae Dong Noh 3 

1 Department of Physics, Korea Advanced Institute of Science and Technology, Daejeon 305-701, Korea 
2 Amari Research Unit, RIKEN Brain Science Institute, 2-1, Hirosawa, Wako, Saitama 351-0198, Japan 
3 Department of Physics, University of Seoul, Seoul 130-743, Korea 
(Dated: February 6, 2008) 

We study the non-equilibrium phase transition in a model for epidemic spreading on scale-free 
networks. The model consists of two particle species A and B, and the coupling between them is 
taken to be asymmetric; A induces B while B suppresses A. This model describes the spreading 
of an epidemic on networks equipped with a reactive immune system. We present analytic results 
on the phase diagram and the critical behavior, which depends on the degree exponent 7 of the 
underlying scale-free networks. Numerical simulation results that support the analytic results are 
also presented. 
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Network concept has been emerging as a useful tool 
for the study of complex interconnected systems 0, 0, Q • 
It allows us to study structure and dynamics of those 
systems in various disciplines in a simple unified manner. 
Many networks in nature share an intriguing property 
of the power-law degree distribution P(k) ~ fc" 7 with 
the degree exponent 7, where P(k) is the probability of 
a node having k links. A network with this property is 
called a scale-free (SF) network pj, examples of which 
include the world-wide web Q|, the Internet and the 
human sexual contact network ®, and so on. 

The broad degree distribution of the SF network brings 
about various interesting phenomena. Particularly we 
are interested in epidemic spreading. Networked systems 
are susceptible to epidemics, which may be computer 
viruses in computer networks or infectious diseases in so- 
cial contact networks of individuals. Robustness against 
the epidemic spreading is one of the key elements for a 
proper function of networks. Researchers have studied 
the characteristic feature of the epidemic spreading in 
the SF network and the efficient immunization strategy 
suitable for the SF network fA H H Hoi fill . 

There are hubs with a large number of links in the SF 
network. This structure makes it vulnerable to an epi- 
demic since it can spread easily through those hubs. It 
is found that an epidemic even with an extremely small 
spreading rate never ceases spreading in the SF networks 
with 7 < 3 0,0. In those networks, the random im- 
munization where all nodes are immunized or vaccinated 
with an equal but finite probability becomes inefficient. 
Instead, the targeted immunization where hubs are im- 
munized preferentially proves to be efficient 0, 0, . 

In the study of the epidemic spreading, one usually 
considers the dynamics of malicious agents on bare or 
immunized networks. On the other hand, the follow- 
ing examples show that the competing spreading dynam- 
ics between malicious agents and immunizing agents are 
also important. Recently, the worm "Code Red" almost 
paralyzed the whole Internet by flooding it with lots of 
useless packets 13] . After the outbreak, there appeared 
the so-called worm-killer worm "Code Green" which was 



meant to seek out and kill the malicious worm 0. Al- 
though it also flooded the network and did more harm 
than good, it hinted a possibility of contagious vaccina- 
tion. The competing dynamics is also observed inside liv- 
ing organisms. When pathogens invade and spread, im- 
mune cells are stimulated. The awaken immune cells then 
replicate themselves and get rid of the pathogens [l5| . 
The phenomenon of cross immunity between competing 
pathogens is another example |12| . If populations are 
exposed to a certain disease, then they become immune 
to the other for certain pair s of diseases, e.g., Hansen's 
disease and tuberculosis |lq . 

In this work, we introduce a model that mimics the 
competing spreading dynamics and study the phase tran- 
sition it displays on SF networks. The model consists of 
two species (A and B) particles, one of which infects 
nodes and the other of which heals the infection. They 
may represent worms (A) and worm-killer worms (B) in 
computer networks or malicious pathogens (A) and im- 
mune cells (B) in living organisms. 

At each time step particles evolve according to the fol- 
lowing dynamic rule: A particle A annihilates sponta- 
neously with the probability pa ■ Or, with the probability 
qA = 1 — Pa j it branches offsprings to its all neighboring 
nodes, a fraction 1 — A (A) of which is the species A (B). 
A particle B annihilates spontaneously with the proba- 
bility pb- Or it branches an offspring B to one of its 
neighboring nodes with the probability qs = 1 — Pb- A 
particle A is removed with the probability /1 upon contact 
with a particle B at the same node. 

Each node may be empty, occupied by a particle A, 
or occupied by a particle B. Such a node can be inter- 
preted as a healthy or susceptible, infected, or immunized 
individual, respectively. Then the process with pa and 
qA corresponds to spontaneous healing and infection, re- 
spectively. The process with A corresponds to activation 
of an immune system. The spontaneous annihilation of 
B particles takes account of the fact that immunization 
may not be permanent. The parameter \i describes the 
efficiency of the immune system in healing the infection. 

Note that the particles A and B have different 
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branching rules. A particle A branches offsprings to 
all neighboring nodes as in the susceptible-infected- 
susceptible (SIS) model Q, while a particle B branches 
an offspring to one of the neighboring nodes as in the 
contact process (CP) ^tJ- This difference can be neg- 
ligible in networks with the narrow degree distribution, 
but it makes a big difference in SF networks 01 ■ We 
adopt the less-reproductive CP dynamics for the particle 
B since it would be costly to activate an immune system 
in real systems. The model was studied in networks with 
a narrow degree distribution ^^|. As we will see, the 
model shows much richer properties in SF networks. 

The model displays three stationary state phases de- 
noted by (0,0), (0,B), and [A,B). They are distin- 
guished with the stationary state density of each particle 
species. Both species are inactive with zero density in 
the (0,0) phase. In the (0,73) phase, the A particles are 
inactive while the B particles are active with a nonzero 
density. Both species are active with nonzero densities in 
the (^4, B) phase. Due to the asymmetric nature of the 
interaction, one does not have a phase in which the A 
particles are active and the B particles are inactive. All 
nodes are healthy in the (0, 0) and (0, B) phases. The 
healthy state requires that the B particles are active in 
the (0, B) phase, which is not the case in the (0, 0) phase. 
In this sense, the (0, 0) phase is the ideal one since the 
remnant B particles would cost system resources. 

We will study the phase transitions and the critical 
behaviors in SF networks by adopting the rate equation 
approach which proves useful for a single species or multi 
species particle systems 0, El ■ It is assumed that the 
particle density at a node is given by a function of its 
degree only. Let a k and b k be the density of A and B 
particles on nodes with degree k, respectively. Then, it 
is straightforward to show that they satisfy the coupled 
rate equations 

a k = -PAak + g,4(l - \)k(l - a k )<d A {k) - ^ika k Q B (k) 
h = -p B b k + q B k{l - b k )e B {k) + Xk(l - b k )Q A {k), 

(1) 

where Q A (k) (@s(fc)) denotes the probability that a node 
with degree k is infected by a particle A (B) from each 
of its neighboring nodes. They are given by @A(k) = 
•£ k ,a kl P(k'\k) and Q B (k) = ^ b k ,P(k'\k)/k' , respec- 
tively. Here P(k'\k) is the conditional probability that a 
node at one end of a link have the degree k' under the 
condition that the other end of the link has the degree 
k. The different form of 9^ and <d B is due to the differ- 
ent spreading dynamics of the two particle species. The 
conditional probability can take care of a degree correla- 
tion |2*oj . In this work, we only consider networks with 
no degree correlation, that is, P(k'\k) = k'P(k')/(k) with 
the deg ree distribution P(k) ~ fc~ 7 and the mean degree 
(k) |2(|. Then, one obtains that 
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FIG. 1: The schematic phase diagram in SF networks with 
7 > 3 in (a) and 7 < 3 in (b). The phase boundary MP 
changes its shape as 7 varies. In (a), it is tangential to MR 
for 7 < 4 (solid line) and non-tangential to MR or MQ for 
7 > 4 (dashed line). In (b), it is parallel to pa axis at M 
for 7 > 5/2 (solid line), and to ps axis near M for 7 < 
5/2 (dashed line). 



They are related to the particle density and will be used 
as the order parameter. 

In the stationary state (a k = b k = 0), the order pa- 
rameter satisfies the coupled self-consistency equations 



q A e A k 2 P(k) 
- {q A Q A + jlQ B )k 



, (3a) 



e B =g{® A ,® B ) = 



(k) 
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(A0A + g B e B )fcP(fc) 
(A6a + q B B )k 



(3b) 



where q A = ^(1 - A), a = A = 3A\ , and q B = 

Note that, without the interaction (A = p, = 0), the self- 
consistency equations are decoupled, each of which has 
been studied separately [HEJ. 

The self-consistency equations allow us to determine 
the phase diagram, whose schematic plot is shown in 
Fig. 2] Let us summarize the result first: For 7 > 
3 (Fig. ma)), we find that the system exhibits the three 
phases (0,0), (0, B), and (A,B). Across the line MR, 
A particles remain inactive and B particles become ac- 
tive. Across the line QM, both particles become active 
simultaneously. Across the line MP, B particles are al- 
ready active and the A particles become active. The 
three phase transition lines merge into the multicritical 
point M. 

On the other hand, the system exhibits the only two 
phases (0,73) and (A, B) for 7 < 3 (Fig. QJb)). Without 
the B species, the A species would always be in the active 
phase for 7 < 3 [?J. The existence of the phase (0,B) 
implies that one can prevent an epidemic from spreading 
even in the SF networks with 7 < 3. The nonexistence of 
the phase (0,0), however, implies that the immunization 
is possible only when the species B is kept to be in the 
active state, which may cost system resources. 

Now we sketch briefly the way the phase diagram 
is obtained, details of which will be presented else- 
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where 21] . It is obvious that Qa = Ob = is a so- 
lution of the self-consistency equation. One can obtain 
the boundary of the phase (0, 0) from the condition for 
the existence of nonzero solutions for or 9^. From 
Eq. we find that /(0,0) = and f (G A , 0) < 1 and 
9 2 /(6a, 0)/dQ 2 A < for all Qa- So a nonzero solution 
for Q A exists when <9/(0, 0)/dO A = q A {k 2 )/{k) > 1 with 
(k 2 ) the second moment of the degree. This gives the 
phase boundary QM at qA = (k) / (k 2 ) or at pa = Va 
with 



Pa 



(l-A)(fc 2 )/(fc) 
l + (l-\)(k 2 )/(k) 



(4) 



Similarly, a 

d g (o 1 o)/de B 



nonzero solution for Q B exists when 
qs > 1, which gives the phase boundary 
MR at q B = 1 or at p B = p c B = 1/2. For 7 < 3, (k 2 ) is 
infinite, p c A = 1, and the phase (0,0) vanishes. 

At the phase boundary MP, ©a = and 9s has a 
nonzero value Q* B satisfying 8 B = g(0,Q B ). We also 
have that 9/(0, Q b )/8Qa — 1 since 0,4 begins to devi- 
ate from zero at that line. These two equations define 
the phase boundary MP. Near the point M, the phase 
boundary is given by 



(-B 



i/mm(i,7-3) 
£ A > 

( 7 -2)/(3--y) 
e A > 



7 > 3, 
7<3, 



(5) 



for small tA—p\— Pa and e B = P B — Pb |21| . 

We also studied the critical behavior of the order pa- 
rameter near the phase transitions by analyzing the prop- 
erty of the functions / and g at small values of Qa and 
Q B . We found that the order parameter shows the power- 
law scaling 
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(6) 



with the 7 dependent critical exponents (3a and (3 B for 
7 ^ 3. The critical exponents have different values at 
different critical lines. The results are summarized in Ta- 
ble[I] The model displays interesting multicritical behav- 
iors at the multicritical point M, which will be discussed 
elsewhere ETA. 

It is noteworthy that the model displays the peculiar 
critical behaviors at 7 = 3. Let us take the degree distri- 
bution as P{k) — ck~ 3 for ko < k with a normalization 
constant c and a degree cutoff kg. With the continuum k 
approximation in Eq. , the functions / and g are given 
by / = ^ln(M^) and g = Y — ^Y 2 ln(I±^) 

with X = q A ®A + A©s and Y = \<d A + q B Q B - The 
logarithmic dependence leads to the following peculiar 
behaviors [2lj: Let = 1 — pa and e B — 1/2 — p B be 
the deviation from M at pa = 1 and p B — 1/2. The 
phase boundary MP is given by 



£A 



(7) 



with a constant d for small 6a and e B . When one ap- 
proaches M from the phase (A, B), the order parameter 
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1 
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e B 


~e- d/eA /e A 
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TABLE I: The exponents Pa and (3b associated with each 
phase transition line. For 7 < 3, the phase (0, 0) denotes the 
line with p A = 1 and pb > Pb- 



exhibits a path dependent critical behavior. Along the 
path with e B = that is tangential to MP, we find that 



Qa 



-2d/e A 



-d/eA 



Along non-tangential paths with finite 6a/\£ B \, we find 
that 



Qa - Oj 



(9) 



In contrast to the power-law scaling at 7 7^ 3, the or- 
der parameter as well as the phase boundary has the 
essential singularity at 7 = 3. The order parameter also 
has the essential singularity near the line with pa = 1 
and p B > p c B (see Table QJ. That is, the transition is 
the infinite order transition. The infinite order transition 
was reported recently in the percolation problem and the 
equilibrium Ising model on growing networks 22, 23] . 

In order to confirm the phase diagram, we performed 
numerical simulations on SF networks generated by us- 
ing the so-called static model j24|. Initially networks are 
filled with the particles, and the density pa and p B of 
each particle species is measured and averaged over at 
least 2000 simulations and 10 network realizations. All 
simulations were performed with A = /i = 0.5. All nu- 
merical data presented below were obtained on the static 
model networks with N = 64 x 10 4 nodes and the mean 
degree (k) = 12. 

We present the numerical data from the SF networks 
with 7 = 3.5 in Figs. [21 (a) and (b). They show the sig- 
nature of the phase transition from (0, 0) to (^4, B) at 
p\ ~ 0.915 and (0, 0) to (0, B) at p% ~ 0.478, respec- 
tively. These numerical data confirms the existence of 
the three phases for 7 > 3. 

At 7 = 2.75, we performed the simulations at p B = 0.3 
and for several values of pa in order to examine the exis- 
tence of the (0, B) phase. The numerical data presented 
in Fig. clearly show that the A species is in the ac- 
tive phase at pa = 0.7 and in the inactive phase at 
Pa = 0.9. Although it is hard to locate the critical point 
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FIG. 2: Density decay at 7 = 3.5. In (a) pA is plotted at 
Pb = 0.9, and in (b) ps is plotted at pa = 0.99. 
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FIG. 3: Density decay on SF networks with 7 = 2.75. pA is 
plotted for several values of pa with ps = 0.3. 

accurately, the numerical data supports the analytic pre- 
diction that the A species can be inactive even for 7 < 3 
in our model. At 7 = 2.75, we also performed the simu- 



lations at pb = 0.9 and at several values of pa close to 
1 in order to examine whether the phase (0, 0) exists or 
not. Numerical data obtained on the networks of sizes 
up to N — 64 x 10 4 nodes indicate that the A species 
is in the active phase at least up to pa — 0.99 and that 
the A species become more active as N increases j2jj ■ So 
we conclude that the phase (0, 0) does not exist in the 
asymptotic N — > 00 limit. 

In summary, we have studied the two-species epidemic 
model on SF networks. The two species A and B are cou- 
pled asymmetrically in that the former induces the latter 
whereas the latter suppresses the former. Our model is 
aimed at describing the spreading dynamics of competing 
malicious pathogens (A particles) and reactive immuniz- 
ing agents (B particles) on complex SF networks. The 
model is shown to have the different phase diagram de- 
pending on whether 7 > 3 or 7 < 3 (see Fig. Our 
results show that one can prevent the epidemic from pre- 
vailing even in SF networks with 7 < 3. The results 
also show that it requires that the immunizing agents 
should be kept in the active phase for 7 < 3, which is 
not necessary for 7 > 3. We have also investigated the 
critical behaviors associated with the phase transitions. 
Especially, when 7 = 3, the phase transitions are infinite 
order transitions with the essential singularity in the or- 
der parameters. 
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